Directed flow of transported and non-transported protons in Au+Au collisions from 

UrQMD model 
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The directed flow of inclusive, transported and non-transported (including produced) protons, as 
well as antiprotons, has been studied in the framework of Ultra-Relativistic Quantum Molecular 
Dynamics approach (UrQMD model) for Au+Au collisions at ^Jsnn =7.7, 11.5, 19.6, 27, 39, 62.4 
and 200 GeV. The rapidity, centrality and energy dependence of directed flow for various proton 
groups are presented. It is found that the integrated directed flow decreases monotonically as a 
function of collision energy for ^/snn =11.5 GeV and beyond. However, the sign-change of directed 
flow of inclusive protons, seen in experimental data as a function of centrality and collision energy, 
can be explained by the competing effect of directed flow between transported and non-transported 
protons. Similarly the difference in directed flow between protons and antiprotons can be explained. 
Our study offers a conventional explanation on the cause of the vi sign-change other than the 
antiflow component of protons alone which is argued to be linked to a phase transition. 

PACS numbers: 24.10.Lx,25.75.Ld,25.75.Nq,24.85.+p 
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I. INTRODUCTION 

Relativistic Heavy Ion collisions, involving two nuclei 
moving close to the speed of light, produce thousands 
of particles. The collective motion of the outgoing par- 
ticles is of special interest because it carries important 
information about the bulk properties of the sytem. A 
handle into the study of such collective motion makes 
use of the azimuthal distributions of detected particles 
in non-central collisions For a quantitative descri 
tion, such distribution is conventionally decomposed 
by a Fourier series 
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where <f> denotes the angle between the particle's az- 
imuthal angle in momentum space and the reaction plane 
angle. The reaction plane is defined by the collision axis 
and the line connecting the centers of two nuclei. The 
various coefficients in this expansion can be calculated 
as 0: 
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The first and the second coefficients are refered to as 
directed (ui) and elliptic flow (V2), respectively, and 
they both play important roles in describing the collec- 
tive expansion in azimuthal space. Elliptic flow is pro- 
duced by the conversion of the initial coordinate-space 
anisotropy into the momentum-space anisotropy, due to 
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the developed large in-plane pressure gradient 0-(5|. Di- 
rected flow, which is the focus of this study, describes the 
"side splash" of particles measured off mid-rapidity @. It 
probes the dynamics of the system in the longitudinal di- 
rection. The shape and magnitude of directed flow in the 
vicinity of midrapidity are of special interest, specially 
for protons, because they are sensitive to the equation of 
state (EOS) and may carry phase transition signal0-Q- 
Directed flow has been measured at the AGS [Io|, IH| , 
at the CERN SPS 0, and at RHIC In par- 

ticular, the directed flow of protons at RHIC has been 
studied in detail by the STAR Collaboration, for Au+Au 
collisions at ^/sWn = 200 GeV [H and at lower ener- 
gies [l8| . It is worthy of notice, that at energies studied 
by STAR, protons and antiprotons exhibit different v\(y) 
slope in close-to-central collisions (5-30% at 200 GeV and 
10-40% at low energies), and at low energies the proton 
vi(y) slope changes its sign from negative as in periph- 
eral collisions to positive as in close-to-central collisions. 
The sign-change of v\(y) slope has been proposed as a 
possible signal of phase transition, it is thus necessary 
to identify other causes that may give rise to the sim- 
ilar sign-change of directed flow in the vicinity of y=0. 
Indeed it is mentioned in [l7j that the observed v\ for 
protons is a convolution of the directed flow of produced 
protons with that of protons transported from beam ra- 
pidity. Since transported protons come from the original 
projectile and target nuclei, they preserve the directed 
flow of spectators. Produced protons, on the other hand, 
tend to flow with the bulk. The flow direction of spec- 
tators and bulk are not necessarily equal; the final v\ 
observed would then be the result of a competition be- 
tween the two different proton sources. 

In this work we study the interplay between the di- 
rected flow of transported protons and that of non- 
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FIG. 1: Particle production in a string-excitation scheme in 
UrQMD. White circles, black circles and gray circles stand 
for initial quarks, produced quarks and produced antiquarks, 
respectively. In case (1), no quark pair is created. In case (2), 
one diquark-antidiquark pair and one quark-antiquark pair 
are spontaneously created in the color flux-tube between the 
initial quark and the initial diquark. In case (3), similarly, 
two quark-antiquark pairs and one diquark-antidiquark pair 
are created, and one meson is produced. 



transported protons, within the framework of the 
Ultra-relativistic Quantum Molecular Dynamics model 
(UrQMD) 0. Note that although vi(y) calculated 
with UrQMD does not reproduce the measurement from 
Au+Au collisions at t/sn~n = 200 GeV, as already 
pointed out in [l7j , the study of the contributions from 
transported and non-transported protons to the overall 
proton v\ can offer insights into the competing dynamics 
which lead to the final observation of the sign of Vi(y) 
and its slope in the vicinity of midrapidity. 



II. TRANSPORTED AND 
NON-TRANSPORTED PARTICLES IN THE 
URQMD MODEL 



second case, one diquark-antidiquark pair and one quark- 
antiquark pair are created. One of the initial quarks 
(white circle) combines with a newly produced diquark 
(black circles) to form the baryon B, the newly pro- 
duced anti-diquark (brown circles) combines with a newly 
produced anti-quark to form the antibaryon C, and the 
newly produced quark together with the initial diquark 
form another baryon D. The third case is similar to the 
second one, and it shows the produced baryons and a me- 
son. At the end baryon A has three initial quarks (3 iq, 
where iq stands for "initial quark" ) , baryon D and H each 
have two initial quarks (2 iq), baryon B has one initial 
quark (1 iq), and baryon G has zero initial quark (0 iq). 
The two anti-baryons (C and G) each have zero initial 
quark (0, iq). The more initial quarks can be found in 
a particle, the more initial information that particle car- 
ries. By tracing the number of initial quarks in a particle, 
we classify particles as produced or transported parti- 
cles. In this article, baryons with three initial quarks are 
regarded as transported baryons, and all other baryons 
that have at least one produced quark are regarded as 
non-transported. A particular subgroup of the latter are 
the produced baryons with zero initial quarks. Produced 
baryons should be similar to anti-baryons in many re- 
spects as both groups are made of produced quarks. 

After the process of fragmentation excitation of 
baryon-strings, baryons produced (mostly unstable 
baryons) are subject to multiple scatterings. When the 
energy of a binary collision is lower than 5 GeV (\/i < 5 
GeV), there will be no strings involved and there are no 
new quark pairs produced. During this process the un- 
stable baryons do decay but the number of constituent 
quarks does not change. For example, in the case of a 
A decay into a proton and a pion, one can still identify 
how many quarks in the daughter proton are originated 
from initial quarks by tracing the quark constituent of 
its parent, A. 



III. ANALYSIS AND RESULTS 



UrQMD is a microscopic transport model based on the 
covariant propagation of all hadrons along classical tra- 
jectories in combination with stochastic binary scatter- 
ings, color string formation and resonance decay. This 
microscopic transport model describes the phenomenol- 
ogy of hadronic interactions at low and intermediate 
energies (s/s < 5 GeV) in terms of interactions be- 
tween known hadrons and their resonances. At higher 
energies (y/s > 5 GeV), the excitation of color strings 
and their subsequent fragmentation into hadrons domi- 
nates the multiple production of particles in the UrQMD 
model 

Fig. [1] schematically sketches the exitation and subse- 
quent fragmentation (y/s > 5 GeV) of a baryon-string 
[19j : In the first case, shown on top of the figure, the 
exitation energy is too low to produce new quark pair(s), 
and there is only one baryon A in the final state. In the 



Fig. [5] shows the rapidity distribution of protons and 
anti-protons for Au+Au collisions at ^/snn = 11.5 GeV 
and 200 GeV, calculated with the UrQMD model. The 
selection on centrality (10%- 70%) and transverse momen- 
tum pt (0.4—1 GeV/c) are chosen to match those used in 
[l7l |. The distribution of transported, non-transported, 
and produced protons are drawn separately. As ex- 
pected, the transported protons peak at beam rapidity 
and the produced ones, at midrapidity. Not surprisingly, 
the distribution of produced protons is similar to that 
of anti-protons (to first order) , and the remaining differ- 
ence comes from their separate production mechanisms 
. For example, for each produced baryon there is an ac- 
companying anti-baryon produced but that statement is 
not necessarily true for anti-baryons. In addition, an- 
tiprotons are subject to annihilation while protons are 
not. The transported (produced) protons accounts for 
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FIG. 2: UrQMD simulation of rapidity distribution of in- 
clusive, transported, non-transported, produced, and anti- 
protons. The simulation is made for 10%- 70% central Au+Au 
collisions collisions at ^/snn = 11.5 GeV (left) and 200 (right) 
GeV. 
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FIG. 3: Rapidity dependence of Vi for protons and anti- 
protons in Au+Au collisions at ^snn = 11.5 GeV for cen- 
trality 10%-70%, from UrQMD calculation. See text for ex- 
planation on individual protons groups. 



12% (2.1%) of inclusive protons at 11.5 GeV, and 0.91% 
(50%) at 200 GeV. The baryon stopping effect is seen at 
both energies. At midrapidity, transported protons dom- 
inate over produced protons at 11.5 GeV, while at 200 
GeV, the opposite is true. Note that the sum of the yield 
of transported and produced protons does not give the 
inclusive proton yield because some of the protons can 
be regarded neither as transported nor produced. Those 
are protons that have 2 or 1 initial quarks, and they ac- 
count for a large fraction of the total number of protons. 
By definition the inclusive proton yield is a sum of trans- 
ported and non-transported proton yields. 

To illustrate the effect of the mixture of initial and pro- 
duced quarks on the shape of v\{y) Figure |3] shows v±(y) 
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FIG. 4: The v\ slope as a function of centrality for protons and 
anti-protons for Au + Au collisions at ^/s^n = 5 11.5 (left) 
and 200 (right) GeV, calculated with the UrQMD model. 



plotted for all proton (antiproton) groups for Au+Au col- 
lisions collisions at ^/snm = 11.5 GeV. It is seen that for 
transported protons (3 iq), the v±(y) slope is positive at 
midrapidity. With the replacement of even only one ini- 
tial quark by a produced quark, as shown by the case 
(2 iq), the v\{jj) shape becomes extremely flat at midra- 
pidity. As more initial quarks are replaced, the Vi(y) 
slope changes sign to negative and continues to decrease 
monotonically with descrcasing number of initial quarks. 
For produced protons (0 iq), the sign of the v\(y) slope 
agrees with that of antiprotons as expected because both 
are produced particles. The v\(y) of inclusive protons 
is a convolution of v\{y) and the relative yields of each 
individual proton group. For clarity, the figures shown 
below do not display the two intermediate proton groups 
(1 iq and 2 iq groups). 

In Fig. @] we present the centrality dependence of v\ (y 1 ) 
slope at midrapidity in Au+Au collisions at y/sJjN = 
11.5 GeV and 200 GeV, for transported, non-transported, 
produced, and inclusive protons. The slope for antipro- 
tons is also shown for comparison. Here the slope is 
extracted from the normalized (y = y/ybeam) rapidity 
distribution, where ybeam is the beam rapidity. Our re- 
sults are obtained by integrating over transverse momen- 
tum (p T ) from 0.4 to 1.0 GeV as done in [l7j|. At 11.5 
GeV, as the overlap between colliding ions range from 
central to peripheral, the vi(y') slope of transported pro- 
tons increases rapidly, reaches a maximum at ^20% cen- 
trality and then decreases. The opposite trend is seen 
for produced protons. The slope of v\{y') extracted 
from non-transported protons is negative in most cen- 
tralities, ranges between positive and negative values, as 
non-transported protons also include intermediate proton 
groups other than produced protons. The v\(y') slope of 
inclusive protons is the result of competing effects from 
directed flow between transported and non-transported 
protons. Because the yield of non-transported protons 
dominate over that of transported protons, in periph- 
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behavior as the one extracted from antiprotons, while the 
Vi(y') slopes for transported protons are all positive. The 
vi(y') slope for inclusive protons is a convoluted effect 
produced by transported and non-transported protons. 
It is positive at 7.7, 11.5,19.6, 27 and 39 GeV, decreases 
monotonically with increasing energy for y / s]v7V =11.5 
GeV and beyond, and becomes negative at 62.4 GeV. The 
sign change of v\{y') slope is also seen in experimental 
data [l8|], although the energy at which the v\{y') slope 
changes sign is lower (11.5 GeV). These findings suggest 
that the apparent "collapse" of v\ , defined by a flat shape 
of the proton vi(y'), or the sign-change of the proton 
vi(y') slope, can be explained by the interplay between 
transported protons and non-transported protons. 
IV. SUMMARY 



FIG. 5: UrQMD calculation of protons' and anti-protons' vi 
slope (dvi/dy ) around mid-rapidity (\y | < 0.5) as a function 
of incident energy. 

eral collisions the sign of the Vi(y') slope agree with that 
of non-transported protons. As the centrality decreases, 
transported proton v\{y') slope increases. This effect 
gradually overcomes the effect of directed flow from non- 
transported protons and change the sign of the vi(y') 
slope to positive. In most central collisions, v\ from all 
groups vanishes due to the symmetrical shape of the col- 
lision. This explains the sign change of v±(y') seen in 
data (3- Similarly, the trend of the v\(y') slope at 200 
GeV can be explained with similar arguments. Note that 
at 200 GeV, the Vi{y') slopes are negative at all central- 
hies. This can be due to the fact that at 200 GeV the 
transported proton yield is too small and the overall ef- 
fect is driven by non-transported protons. For the same 
reason, a difference between proton and antiproton v\ is 
observed as seen in the inset, similar to what is shown in 

The energy dependence of directed flow for protons and 
anti-protons is presented in Fig. [5] Based on UrQMD 
calculation, the slopes of v\(y') extracted from produced 
protons are negative and they follow more or less a similar 



Directed flow of inclusive, transported protons and 
non-transported protons (including produced protons), 
as well as antiprotons, have been studied with the 
UrQMD model at RHIC energies. It is found that the in- 
tegrated directed flow decreases monotonically as a func- 
tion of collision energy for ^/snn =11.5 GeV and be- 
yond. However, the sign-change of directed flow of inclu- 
sive protons, as a function of centrality and collision en- 
ergy, can be explained by the competing effect of directed 
flow between transported and non-transported protons. 
Similarly, the difference in directed flow between protons 
and antiprotons can be explained. Our study identified 
alternative cause of the vi sign-change other than the 
anti-flow component of protons alone which is argued to 
be linked to a phase transition. 
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